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M.I.T. FLOATING POINT ARITHMETIC PACKAGE 



Program Library Write-up 



DECUS No. 84 



The Floating Package is a group of arithmetic subroutines in which 
numbers are represented in the form f x 2 e . f is a one's complement 
18-bit fraction with the binary point between bits and 1. e is a 
one's complement 18-bit integer exponent of 2. The largest magnitude 
numbers that can be represented are ^ Q 39,000. 

A number is normalized when x ^ /f / < 1. All the floating routines, 
except the two floating unnormalized adds, return a normalized 
answer. The fraction appears in the ac, the exponent in the io. 
Description of routines: 

Floating Add - jda fad 

One argument should appear in the ac-io. The other argument should 
have the addresses, direct or indirect, of the fraction and exponent 
in the two registers following the jda fad . 



lea fl 
lio el 
jda fad 

f2 

e2 

dac 



/load first argument 

/call floating add 
/address of second fraction 
/address of exponent for second 

fraction 
/control returns to here with 
/normalized answer in ac-io 



Floating Multiply - jda fmp 

lac fl 
lio el 
jda fmp 

f2 

e2 
dac 



/Load multiplicand. 

/Call floating multiply 
/address of fraction of multiplier 
/address of exponent of multiplier 
/Control returns to here with normal - 
/ized answer in ac-io. 



Floating Divide - jda fdv 

lac fl 
lio el 
jda fdv 

f2 

e2 
dac 



/Load first argument. 

/Call divide. 

/address for divisor, hit will occur 

/if f2=0. 
/Control returns here with normalized 
/answer in ac-io. 



-2- 



Floating Square Root - jda fsq 

Execution time '^385 /^sec. 

lac f 

lio e 
jda fsq 

dac 



Floating Log, base 2 - jda log 

lac f 
lio e 
jda log 
dac 



Floating Reciprocal - jda rep 

lac f 
lio e 
jda rep 

dac 



Floating Input - jda fip 

jda fip 
jsp 



dac 



Legal characters for fip 



/Load argument; argument must be 
normalized 

/Call square root; hit will occur 
/unless fZO. 

/Control returns here with normalized 
/answer in ac-io. 



/Load argument. 

/Call log; hit will occur unless f > 0, 
/Control returns here with normalized 
/answer in ac-io. 



/Load argument. 

/Call reciprocal; hit will occur in 
/fdh if f=0. 

/Control returns here with normalized 
/answer in ac-io. 



/Call input; ac-io don't matter. 
/This instruction is repeatedly 
/executed ( xct ) in order to get the 
/input characters. The jsp (or jda ) 
/could call a typewriter or reader 
/listen loop subroutine which should 
/return the input characters in the 
/low bits of the io. 

/Control returns here with the answer 
/in the ac-io after the first illegal 
/character. 



x resets routine and starts forming a new number. Spaces 
and code deleted characters are ignored. Legal characters 
are: ., e, 0-9, -, x, and space. The illegal character 
that terminated the number is in register fip . 



-3- 



Input examples: 

6.9el 

690 e-1234 

-6.9 e 17 

Floating Output - jda fop 

lac f 

lio e 

j da fop 
tyo 

dac 



/Load argument. 

/Call output. 

/an executed instruction ( xct ) for 

/output 

/Control returns here with normalized 

/floating point input quantity. 



The routine generates 
put instruction could 

The output format is 

Floating Unnormalized 

lac fl 

lio el 

jda fua 

f2 

e2 



parity for each character, so the executed out- 
be a ppa or a call to an output subroutine. 

71000 e2, 5 significant figures. 

Add - jda fua 

/Load first argument. 

/Call unnormalized add. 
/addresses of second argument 



The subroutine returns with a 35 bit number in the ac-io with binary 
point after the bit number equal to the larger exponent of the two 
arguments. If the addition produces an overflow, the larger exponent 
is incremented by 1. In any case, the larger exponent, perhaps 
incremented, appears in fac+1. Examples for subroutine: 



lac (200000 A 
lio (0 
jda fua 

(0 

(17. 

At return ac,io equals 0,400000. 

lac (0 
lio (16. 
j da fua 

(200000 

(0 



/zero with exponent to cause the number 
/to be fixed. 



/0 



/- 



At return ac,io equals 1,0 



-4- 

Floating Unnormalized Add and Round - jda fur 

lac fl /Load argument 

lio el 

jda fur /Call unnormalized add and round 

f2 /addresses of second argument 

e2 

dac /Control returns here with fraction 

/in ac and exponent in io. 

This routine is the same as fad except that the answer is not normal- 
ized. The larger exponent returns in the io, unless overflow occured. 
Then the larger exponent +1 returns in the io. Example: 

lac (300000 /3 

lio (2 

jda fur /Call subroutine. 

(0 /zero with exponent to 

(17. /cause the answer to be fixed. 

At return, ac,io equals 3,17. 

Floating Exponentiation - jda f2x 

This subroutine calculates 2 X . Execution time '2" 1.3 m sec. 

lac f /Load argument, 

lio e 

jda f2x /Call subroutine. 

dac /Control returns to here with normalized 

/answer in ac-io. 



/floating pack with trig func, arc tan, and nat log 



fac, 










fmp, 







/floating multiply 




dap 


fml 






dap 


fm2 






idx 


fm2 






dao 


frx 






idx 


frx 






swp 






fm2, 


add 


i • 






dac 


fac 1 




fml, 


lac 


i • 






mul 


fmp 






dac 


fac 






jmp 


f nm 





fnm, cla /normalize. For internal use only 



flp, 



fdf, 
fze, 



fnr, 



frn, 



frx, 



cla 


/normalize. 


dap fnr 




add fac 




sza i 




jmp fze 




sub (200000 


sma 




jmp fnr 




idx fnr 




lac fac 




scl 1 




dac fac 




jmp flp 




law 17. 




dap fnr 




lac fac 




scl 9s 




scl 8s 




add (0 




sza 




jmp fdf 




law i . 




add fac 


1 


dac fac 


1 


lac fac 


/round • 


sar 9s 




sar 8s 




clo 




scl 1 




add fac 




lio fac 


1 


szo 1 




jmp • 




rar 1 




swp 




idx fac 


1 


swp 




jmp frx 





fdv, /floating divide 

dap fdl 

dap fd2 

ldx f d2 

dap f rx 

ldx f rx 

swp 
fd2, sub 1 . 

dac fac 1 

lac fdv 

ell 

spa 

llo (-0 
fdl, dlv 1 . 

Jmp fdo 

dac fac 

cla 

spl 

cma 

swp 

dlv 1 fdl 

hit /impossible 

scr 9s 

scr 8s /or mul (l 

Jrap f nm 
fdo, add (0 

sza 1 

hit /zero divisor 

scr 1 

dac fdv 

ldx fac 1 

lac fdv 

Jrap fdl 



J 



fad, 



fzn, 



fnz, 



fan, 








dap 


fal 


dap 


fa2 


idx 


fa2 


dap 


f rx 


idx 


frx 


law 


fnm 


dap 


fux 


cla 




add 


i fal 


sza 




jmp 


fnz 


lac 


fad 


dac 


fac 


dlo 


fac 1 


cli 




spa 




lio 


(-0 


jmp 


fnm 


cla 




add 


fad 


sza 




jmp 


fah 


lac 


1 fal 


lio 


i fa2 


jmp 


fzn 


dac 


fua 


jmp 


fun 



/floating add 



/ 



fua, 





/unno 




dap 


fal 




dap 


fa2 




idx 


fa2 




dap 


fux 




idx 


fux 


fun, 


swp 






dac 


fac 1 


fa2, 


sub 
sma 


i . 




jrap 


fsr 




dac 


fde 




lac 


i fa2 




dac 


fac 1 




lac 


fua 




dac 


fa4 


fal, 


lac 


i • 




dac 


fua 


fba, 


cli 
spa 






lio 


(-0 




dio 


fa4 1 




scr 


1 




and 


(377777 




dac 


fmp 




swp 






sar 


A 
-A- 




and 


(311111 




dac 


fad 




lac 


fde 




scr 


3s 




add 


(5 




spq 






jmp 


f ou 




add 


.+1 




da p 


fsh 




clc 






scl 


3s 




add 


(fsp 7 




dap 


fxq 




lac 


fa4 




cli 






spa 






lio 


(-0 


fsh, 


jmp 


• 




repeat 4, scr 8s 


fxq, 


xct 


• 



/unnormalized fleeting add (fix) 



swp 

sar 1 

and (377777 

add fad 

dac fa3 1 

spa 





idx 


fmp 






swp 








oma 








sub 


f mp 






cma 








dac 


fa3 






sma 








jmp 


ful 






law 


i 1 






sub 


fa3 1 






cma 








dac 


fa3 1 






sas 


(400000 






jmp 


ful 






law 


i 1 






sub 


fa3 






cma 








dac 


fa3 




ful, 


lac 


fa3 






lio 


fa3 1 






ril 


1 






rcl 


1 






scr 


2s 






scl 


2s 






dac 


fac 






xor 


fa4 






dac 








lac 


fac 






xov 


fua 






and 


fa3 






sma 








jmp 


fok 






idx 


fac 1 






lac 


fac 






scr 


1 






xor 


(400000 






dac 


fac 




fok, 


lac 


fac 






dlo 


fa4 1 




fux, 


jmp 


• 




fsr, 


cma 








dac 


fde 






lac 


i fal 






dac 


fa4 






lac 


fua 






jmp 


fba 




f ou, 


lac 


fua 






lio 


fa4 1 






jmp 


fok-1 












fc 


fs3=o 








f SD, 


repeat 8., scr 


8s-fs8 


fde, 










fa4, 

fs8=fs8 fs8 1 



fur, /unnormalized rounded add 

dap ful 

dap fu2 

idx f u2 

dap frx 

idx frx 

lac fur 

jda fua 
ful, i . 

fu2, i . 

jmp frn 

/square root jda fsq Inputs must be normalized (or 0) 
fsq, 

dap zlv 

law i 3 dac zlv 1 

spi 

cma 

rcr 1 

dio zlv 2 

x or zlv 2 

swp 

spi sma-skp 

idx zlv 2 

lac fsq 

spa 

hit 

sza i 

jmp zlv-1 

spi 

sar 1 

dac fsq 

sar 2s 

add zlv 3 

jmp . 11 

lac fsq 

cli 

scr 2s 

dis zlv 4 

hit 

add zlv 4 

cli 

rcr 1 

dac zlv 4 

isp zlv 1 

jmp • -12 

lac zlv 4 

scl 1 

sza i 

sub (400000 

lio zlv 2 
z Iv , jmp 



66314 



/log, base 2. requires normalized arg 
log, 



lgo, 










dap lgo 






dio lgo 


6 




lac log 






spq 






hit 






sar 1 






add lgo 


1 




dac lg 






la c 1 og 






sub lgo 


2 




cli 






spa 






lio (-0 






div lg 






hit 




/not norm. 


dac lg 




mul lg 


mul lgo 


3 




sar 4s 






add lgo 


4 




mul lg 






sub lgo 


5 




lio (1 






jda fad 






lgo 6 




(17. 


jmp . 






132405 
270517 




265012 373621 
100002 



/floating 


input! jda 


fip, input inst 


fip, 











dap 


owt 


dap wat idx ow 


ini, 


dzm 


hoi 






dzm 


hoi 1 






dzm 


zll 






dzm 


dig 






lio 


wat 3 


/spa 




dio 


ens 






dio 


enn 






dio 


6fg 




wat, 


xct 


• 






dio 


fip 






rcr 


7s 






spa 








Jmp 


wat 






sar 


2s 






sar 


9s 






lio 


fns 2 


/sma 




sad 


(chara< 


2 rx 




Jmp 


ini 






sza 


i 






Jmp 


wat 






sas 


(chara< 


2 r- 




Jmp 


5fg 






dio 


ens 






Jmp 


wat 




5fs, 


sas 


(chara< 


2 r. 




Jmp 


. 3 






dio 


enn 






Jmp 


wat 






sas 


(charai 


2 re 




Jmp 


cnm 






dio 


6fg 






lio 


wat 3 






dio 


ens 






Jmp 


wat 




cnm, 


sad 
cla 


har 1 






sub 


wah 






sma 








Jmp 


fns 






add 


wah 




ens , 


spa 
cma 


• 






dac 


zl2 






cla 






6fg, 


spa 


• 






Jmp 


cxp 




enn, 


spa 


• 






idx 


dig 






lac 


hoi 






lio 


hoi 1 






jda 


f mp 






wan 








har 


5 






Jda 


fad 






zl2 








har 


5 






dac 


hoi 






dio 


hoi 1 






Jmp 


wat 





cxp, lac zll 

sal 2s 

add zll 

sal Is 

add zl2 

dac zll 

jmp wat 
fns , lac zll 

sub dig 

sma 

llo . 1 

spa 

cma 

dac zll 

dlo cpr 
cpr, spa .-• 

jsp Inv 

law har 

dap tnp 

law har 1 

dap tnp 1 
rpt, lac zll 

sza i 

jmp ard 4 

scr 1 

dac zll 

spl i 

jmp ard 

lac hoi 

lio hoi 1 

jda fmp 

tnp 1 

tnp 1 1 

dac hoi 

dio hoi 1 
ard, Idx tnp 1 

dap tnp 

idx tnp 1 

jmp rpt 

xct cpr 

jsp inv 

lac hoi 

lio hoi 1 
owt, jmp . 

inv, dap qzl 

lac hoi 

sza i 

jmp owt-1 

lio hoi 1 

jda rep 

dac hoi 

dio hoi 1 
qzl, jmp .-. 

wah, 10. 

zll, zl2, 



/floating output 

fop, 

dap urp 

dap xit 

idx xlt 

law har 30, 

dap tnp 

law har 31. 

dap tnp 1 

cla 

dap ubm 

lac fop 

dlo 151 

lio 2hn 

sma 

jmp urp 

cma 

dac fop 

lio (charac r- 
urp, xct •-. 

lac f 00 

lio 151 

Ma fmp 

(5 

har 5 

dzm loh 

sza i spi-skp 

jmp mzd 

jda rep 

dac fop 

rcr 9s 

rcr 9s 

sub (1 

dac 151 

law spa-skp sma-skp 

dap ubm 
mzd, dzm dig 

lac fop 

lio 151 

da c Is t 

dio 1st 1 

jda fdv 
tnp, har 30. 

har 31. 

dac fop 

dlo 151 

jda fmp 

(5 

har 5 spi 

jmp tuf 
idx dig 
jmp mzd 1 



"GUI 



lio 1st 1 



2hn, 



ubm, 



drp, 



xi i, , 
har, 

17. 

97656. 

99565. 

116246. 

100343. 

123712. 

hoi, 
1st, 
cl, 
dig, 



dac 
dio 
lac 
sal 
add 
dac 
law 
add 
dap 
dap 
idx 
sas 
jmp 
skp 
jmp 
lac 
lio 
jda 
fop 

151 
dac 
dio 
lio 
xct 
lac 
lio 
jda 
tn5 
bar 

fo a 

har 
jda 
xct 
lio 
xct 
lio 
xct 
lac 
cma 
xct 
cma 
jda 
xct 
jmp 
20 o 



fop 
151 
loh 
200 

dig 
loh 
i 2 

tnp 
tnp 
tnp : 

tnp : 

tob 

mzd 

• *" . 
drp 
cl 

cl 1 
fdv 



fop 

151 

(charac r, 

urp 

fop 

151 

f mp 



fur 

5 

dpt 
urp 
2hn 

urp 

(265 

urp 
loh 

ubm 

dpt 

urp 



20 



100. 



27. 

213. 

1701. 

13607. 

108853 

na.r— x 





314632 





72759. 54. 

75632. 426. 

103097. 3402. 

76818. 27214. 

-i-uu, ^ 





-3 

151, 



17. 

80778. 

872 83. 

810Q3. 
90042 . 



10000. 

107. 
851. 
6804. 
54427. 



tn5, 



12500. 



*, 



/decimal integer print of ac. jda dpt followed by output instr, 

dpt, 

dap dpo 

dap dpx 

dzm ddv 

idx dpx 

lio (charac r- 
dlp, lac dpt 

spa 

xct dpo 

spa 

cma 

dac dpt 
dll, dac dpr 

mul (1 

div (10. 
dpr, 

sas ddv 

jmp dll 

sni /note sni 

lio (charac rO 
dpo, xct • 

lac dpr 

dac ddv 

sas dpt 

jmp dip 
dpx, jmp . 

ddv, 

/parity for low 6 io bits, saves ac 
Pty, 

dap ytp 

law i 770 

rcr 6s 

lio (252002 

rcr 9s 

dap . 1 

rir 

spi 

and pty 2 

rcl 6s 

rcl 9s 

lac pty 
ytp, jmp . 

/reciprocal routine 
rep, 

dap per 

dio per 1 

cli 

law 1200 

rcl 9s 

jda f dv 

rep 

per 1 
per, jmp . 

per 1, 



/antilog, 
f2x, 



fci, 
fmt, 



(jda rep 

fck 



f2x 

fua 



fmp 

17. 
fac 1 



/power too big 
(200000 
fac 



fzO, 



sni i 
jmp fci 

fdu, lac fac 
lio fmp 

fck, •-. . 

f xx , jmp 



ftb 16. 



ftb, 265017 230160 213454 205526 202633 
201312 200544 200262 200131 200054 200026 
200013 200006 200003 200001 200001 



A to 


the y power 


pow, 







dap pol 




dap po2 




idx po2 




dap pox 




idx pox 




lac pow 




jda log 




jda fmp 


pol, 


1 • 


Po2, 


1 . 




jda f2x 


pox, 


jmp . 



/floating normal angle, for Internal use only 

fna, 

dap rtl 
elf 6 
dio Fed 
lac fna 
spa 

jmp inc 
lac fna 
Ho fed 
jda fad 
[-311040 

(3 

sma 

dac fna 
sma 

dio fed 
sma 

jmp . -12 
lac fna 
lio fed 
jda fad 
[-226630 

(3 

sma 
jmp 4q 
jda fad 
(311040 

sma 
jmp 3q 
.Ida fad 
[311040 

(1 

sma 

jmp 2q 

lac fna 

lio fed 
rtl, jmp . 
inc, lac fna 

jda fad 

[311040 

(3 

dac fna 
dio fed 
jmp fna 4 
4q, dac fna 
dio fed 
stf 6 
lac fna 

lio fed 
2q, cma 

Ida fad 
T3H040 

dac fna 
dio fed 
jmp rtl-2 
3q, dac fna 

dio fed 



stf 6 
jmp rtl-2 

/floating sin-cos 

fcs, 

dap rt2 
lac fcs 
jda fad 
1311040 

(l 

jda fna 

dac fsn 
dio fex 
jmp .+5 
fsn, 

dap rt2 
lac fsn 
jda fna 
dac fsn 
dio Fex 
jda fmp 
fsn 
fex 

dac Fxs 
dio Fes 
jda fmp 
fsn 
fex 
jda fmp 

(252533 
(-2 

dac Fdd 
dio Fed 
cma 

jda fad 
fsn 
fex 

dac fcs 
dio Fee 
lac fdd 
lio fed 
jda fmp 
fxs 
fes 
jda fmp 

{314637 

dac fdd 
dio fed 
jda fad 
fcs 
fee 
dac fcs 

H 1 r\ ■Pan 

VA X V J. ^> \_ 

lac fdd 
lio fed 
jda fmp 
(303034 
(-5 



rt2, 



fsc, 



fco, 



rt3, 



jda 


fad 


fcs 




fee 




sza 


i 


lio 


fex 


sza 


I 


lac 


fsn 


szf 


6 


cma 




jmp 


• 


secant - co: 







dap 


rt3 


lac 


fsc 


jda 


fad 


(311040 


(1 




dac 


fco 


jmp 


.+3 







dap 


rt3 


lac 


fco 


jda 


fsn 


sza 


1 


jmp 


.+3 


jda 


rep 


jmp 


• 


lio 


(377777 


jmp 


rt3 



lal 



// 



/floating tan cot 

ftn, 



rt4, 



fct, 



rt5, 








dap 


rt4 


dio 


Fte 


lac 


ftn 


jda 


fcs 


dac 


fsc 


sza 


i 


jmp 


rt4 1 


dio 


fco 


lac 


ftn 


lio 


fte 


jda 


fsn 


jda 


fdv 


fsc 




fco 




jmp 


• 


lio 


(377777 


lai 




jmp 


rt4 







dap 


rt5 


dio 


fte 


lac 


fct 


jda 


fsn 


dac 


fsc 


sza 


i 


jmp 


rt5 1 


dio 


fco 


lac 


fct 


lio 


fte 


jda 


fcs 


jda 


fdv 


fsc 




fco 




jmp 


• 


lio 


(377777 


lai 




jmp 


rt5 



fin, /floating natural log 



lrp, 



rt6, 








dap 


rt6 


dlo 


eln 


lac 


fin 


spq 




hit 




lac 


(200000 
lnc 


dac 


lac 
dac 


I 1 

Tee 


lac 


fin 


Jda 


fad 


lnc 




Ice 




dac 


fin 


dlo 


eln 


Jda 


fad 


[-200000 


(2 




jda 


fdv 


fin 




eln 




dac 


fin 


dio 


eln 


dac 


Tna 


dio 


Tne 


jda 


fmp 


fin 




eln 




dac 


fls 


dlo 


els 


lac 


lnc 


lio 


Ice 


Jda 


fad 


[200000 


(2 




dac 


lnc 


dio 


Ice 


lac 


fin 


lio 


eln 


jda 


fmp 


fls 




els 




Jda 


fdv 


lnc 




Ice 




jda 


fad 


lna 




lne 




sad 


lna 


Jmp 


rt6 1 


dac 


lna 


dio 


lne 


jmp 


• 


r 


fmp 
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